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The rich physics exhibited by random optical wave fields permitted Hanbury Brown and 
Twiss to unveil fundamental aspects of light. Furthermore, it has been recognized that optical 
vortices are ubiquitous in random light and that the phase distribution around these optical 
singularities imprints a spectrum of orbital angular momentum onto a light field. We demon¬ 
strate that random fluctuations of intensity give rise to the formation of correlations in the 
orbital angular momentum components and angular positions of pseudothermal light. The 
presence of these correlations is manifested through distinct interference structures in the or¬ 
bital angular momentum-mode distribution of random light. These novel forms of interference 
correspond to the azimuthal analog of the Hanbury Brown and Twiss effect. This family of 
effects can be of fundamental importance in applications where entanglement is not required 
and where correlations in angular position and orbital angular momentum suffice. We also 
suggest that the azimuthal Hanbury Brown and Twiss effect can be useful in the exploration 
of novel phenomena in other branches of physics and in astrophysics. 

Introduction 

In 1956, Hanbury Brown and Twiss (HBT) revolutionized optical physics with the observation of a new form 
of interference produced by correlations of the intensity fluctuations of light from a chaotic source. Their 
stellar interferometer collected light produced by independent sources on the disc of a star and detected at 
two different locations on Earth . The observation of a second-order interference effect in this configuration 
was intriguing because at that time it appeared that classical and quantum theories of light offered different 
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predictions [^. Ever since, this effect has motivated extensive studies of higher-order classical correlations and 
their quantum counterparts in optics, as well as in condensed matter and particle physics . Fundamental 
bounds have been established for the degree of correlation for a wide variety of degrees of freedom, such 
as in polarization, time, frequency, position, transverse momentum, angular position and orbital angular 
momentum (0AM) (7j|^. 

The random nature of light is an essential element of the HBT effect. Moreover, the random proper¬ 
ties of light have been investigated and applied in a wide variety of other contexts. For example, speckled 
light, intimately related to pseudothermal light, has played a fundamental role in the development of optical 
physics, imaging science, and nanophotonics. In addition, the study of fundamental processes such as trans¬ 
port phenomena, localization of light, optical vortices, and optical correlations has led to the development 


of novel physics produced as a consequence of the chaotic properties of light 10 -15 . These results have 
motivated interest in the design of random lasers and of disordered structures that scatter light in random 
directions, which serve as sources of pseudothermal light p6] . 

As identified by Berry, optical vortices produced by the interference of random waves are intrinsic el¬ 


ements in chaotic light 17 18 . Interest in this field has exploded since the recognition of a special class 
of vortices that carry 0AM, characterized by an azimuthal phase dependence of the form where t is 


the 0AM mode number and (j) is the azimuthal angle 19 . The azimuthal properties of light, described by 
the conjugate variables of angular position and 0AM, have shown potential for technological applications in 


information science, remote sensing, imaging, and metrology 20 . In astrophysics, recent theoretical studies 
have predicted that rotating black holes can imprint an 0AM spectrum onto light. The measurement of 


this spectrum could lead to an experimental demonstration of the existence of rotating black holes 21 


In addition, the optical vortex coronagraph has allowed the observation of dim exoplanets by canceling a 


diffraction-limited image of a star 22 . More recently, it has been proposed to use rotational Doppler shifts 


for astronomy 23 


Here, we show that random fluctuations give rise to the formation of intensity correlations among the 
0AM components and among the angular positions of pseudothermal light. Furthermore, we show that the 
presence of these correlations leads to a variety of complex interference structures that correspond to the 
azimuthal analog of the HBT effect. In the original HBT experiment, two detectors were used at different 
locations to gain information about the physical size of a distant incoherent source. In our experiment, we 
use two detectors to measure intensity correlations between two 0AM components of an incoherent source 
with controllable spatial and temporal coherence. We show that such correlations unveil the azimuthal 
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structure of the source, which is shaped in the form of double angular slits in our realization. We study the 
far-held pattern by projecting it onto various 0AM modes, and measure hrst- and second-order interference 
patterns of this structure. We identify two key signatures of the azimuthal HBT effect. The hrst is that 
HBT interference can show features in the 0AM mode distribution at both the frequency and at twice the 
frequency of the hrst-order coherence produced by coherent light. The second consists of a shift of the 
interference structure when plotted as a function of 0AM. We hnd that each of these ehects depends on the 
strength of the huctuations of the pseudothermal light. We also study the nature of the correlations between 
different 0AM components and between different angular positions of pseudothermal light, and we find that 
these depend on the strength of the fluctuations as well. These effects correspond to the classical counterpart 
of azimuthal Einstein-Podolsky-Rosen (EPR) correlations [^, and throughout this article, we highlight the 
similarities and differences between thermal and quantum correlations as manifested in the azimuthal degree 
of freedom. 


Results 

Origin of HBT interference in the OAM domain 

As in the original HBT experiment, we collect light from two portions of a random field. This is carried out 
through the use of two angular slits. We represent the optical field after the slits as 

= £{r)<^{r,(j))[A{(j)) + A{(j) - (j)o)]. (1) 

Here, £(r) represents the coherent optical field produced by a laser, $(r, 0) is a particular realization of a 
random phase screen, and A((/)) describes the transmission function of the angular slits. A{cj)) is centered at 
0 radians, and, therefore, A{(p— ((iq) is centered at (j>o- We next consider the projection of the optical field of 
Eq. (1) onto a set of OAM modes. The result of such a measurement is described by the quantity Upg defined 
as J r dr d(j) {2 tt where i?; (r) is a radial mode function with radial index p and i is 

the OAM index. Consequently, the measured intensity for each OAM projection Ig is equal to 
The average of the intensity over an ensemble of different realizations of the fluctuating field is then given 

by 

ih) = J rdr|£-(r)|2{2 + e-*^‘^«(d>*(r,O)d>(r,0o))+e*^‘^«(d>*(r,</>o)$(r,O))}, (2) 

where a is the width of the slits, and the ensemble average is denoted by (...). It is evident that the angular 
double slit gives rise to Young’s (first-order) interference in the OAM-mode distribution of the optical field 
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and that this interference is dependent on the angular separation of the two slits, (j)Q. Furthermore, the 
visibility of the interference pattern is determined by the terms ($*(r, 0)$(r, (/)o)) and ($*(r, (/)o)$(r, 0)), 
which quantify the field correlation between two different angular positions. These terms are sensitive to the 
phase difference of the field at two points. Consequently, the interference visibility becomes smaller as the 
degree as spatial coherence is reduced. 

In direct analogy to the HBT experiment, in which two detectors measure the transverse momentum 
(far-held) distribution of a random held emitted from two locations of a star, we measure the correlation 
between two 0AM components of light emitted from a random source shaped as two angular slits. Similar to 
linear position and linear momentum, angular position and 0AM are conjugate variables and form a Fourier 
pair. Thus, we consider the second-order coherence function G\l= {hilt 2 )^ which is the key quantity that 
describes the azimuthal HBT effect. This quantity is a measure of the intensity correlations between the 
components of the the held with 0AM values and ^ 2 ■ 

We consider a special case in which we measure the the second-order correlation at symmetrically dis¬ 
placed 0AM values of ^ and —1. In the context of the original experiment of HBT, this situation would 

involve measuring the receiving apertures by equal amounts in opposite directions. To analyze this situation, 

f2l 

we need to determine the second-order coherence function G\ = {Iil-i). We hnd that this quantity can 
be expressed (see Supplementary Materials) as 

ihl-i) = Qo + Gi + Gil- (3) 

The intensity correlation function thus consists of three contributions. The hrst is a constant term denoted 
by ^0 whose form is shown in the Supplementary Materials. The second term, Gt^ describes an interference 
pattern that oscillates in I at the same frequency as {Ig) and is given by 

Gi = ^ J ridrir 2 dr 2 \S{ri)\^ l^(^2)|^ (e“*'^'^‘’{($*(ri,0)$(ri,(()o)) -f ($*(r2, 0o)$(?'2,0))} -f c.c.). 

(4) 

The last term, shows an interference pattern that oscillates in the 0AM value ^ with twice the frequency 
of {Ig) and it is given by 

G 2 t = ^ J ridrir 2 dr 2 |■^^(?'l)|^ \£{r 2 )\^ (e“^*^‘^“{($*(ri, 0)$(ri, (/)o)4‘*(?'2, <()o)^’(?'2,0))} -f c.c.). 

(5) 

We see that the contribution Gg depends on a phase-sensitive term (<I)*(r, 0)<I>(r, (/iq)) that decreases in 
magnitude with increasing randomness induced by field fluctuations. The visibility of this contribution to 
the interference pattern thus decreases with increasing field fluctuations. However the contribution G 2 (. is 
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proportional to a positive-definite quantity (|$(r, 0)|^ |<i)(r, 0o)|^) that survives even in the presence of the 
fluctuations in the chaotic field. 



Figure 1: Experimental setup for the study of the azimuthal HBT effect. (A) The 532 nm output of a solid laser is directed 
onto a digital micro-mirror device (DMD), where a random transverse phase structure is impressed onto the beam. A 4/ optical 
system consisting of two lenses with different focal lengths (figure not to scale) and a pinhole is used to isolate the first diffraction 
order from the DMD, which is a pseudothermal beam of light. This beam is then passed through a beam splitter (BS) to create 
two identical copies. Each copy is sent to a separate spatial light modulator (SLM) onto which a computer-generated hologram 
is encoded. (B) For the HBT measurements, a pair of angular slits is encoded onto the SLMs. In addition, forked holograms 
corresponding to 0AM values are encoded onto the same holograms to project out controllable 0AM components. For our 
measurements of the 0AM and angular-position correlation functions, we do not use the double slit but simply project onto 
0AM values or angular wedges, respectively. (C) Intensity distribution of a generated pseudothermal beam of light. 


Experimental demonstration of azimuthal HBT interference 

Our experimental setup is depicted in Fig. 1 (A and B). We use a solid state laser working at 532 nm 
along with a digital micro-mirror device (DMD) and a 4/ optical system containing two lenses and a spatial 
filter in the Fourier plane to isolate one order of diffraction from the DMD. We first impress a sequence 
(at a 1.4-kHz writing rate) of random transverse structures having Kolmogorov statistics onto the beam to 
simulate thermal light [^[^. For details, see Materials and Methods. This procedure modifies the spatial 
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and temporal coherence of the beam in a fashion similar to the modification induced by a rotating ground 
glass plate (see the intensity distribution of the beam in Fig. 1 C), which is often used to produce light 
with thermal statistics. We quantify the spatial coherence of the beam by means of the Fried coherence 


length ro 27 . The strength of spatial phase variations within the beam increases as ro decreases. By virtue 
of ergodicity, iterating through an ensemble of such holograms results in random phase fluctuations in time 
characterized by the parameter rg. The structured beam is then split into two parts at a beam splitter, and 
each is imaged onto a spatial light modulators (SLM). On each SLM, a pair of angular slits and a forked 
diffraction grating are encoded (see Fig. IB). The first diffraction order of the SLM is collected by a single¬ 
mode optical fiber (SMF), measured by avalanche photodiodes (APDs), and their degree of correlation is then 
computed. The time window for determining coincidence events is set to 42 ns, and the total accumulation 
time is set to 15 s. 

We begin with the measurement of first-order (Young’s) interference in the 0AM domain, which can be 
observed in the OAM-mode distribution of light measured by either of the two detectors. For each value of £, 
we impress several hundred random phase screens onto the DMD, all characterized by the same value of rg, 
and we then calculate the correlation of the intensity. We repeat the experiment for all I in the range ^ = —15 
to £ = -1-15. We perform this task by encoding holograms onto the SLMs in which the two angular apertures 
are multiplied by different forked diffraction gratings (see Fig. IB). The OAM-mode distributions of the 
field as given by (/^) are shown in Fig. 2 (A to D). Fig. 2A shows the interference obtained when spatially 
coherent light is used, and Fig. 2 (B to D) shows the interference for different regimes of pseudothermal light, 
as characterized by successively decreasing values of rg. The visibility is seen to decrease with the decrease 
of the spatial coherence of the source. 

We next study second-order coherence. Our experimental results for the second-order coherence function 
AG) defined as Qi + Qu-, are shown in Figs. 2 (E to H). For a coherent beam (Fig. 2E), Qi, is the dominant 
contribution to AG) We reach this conclusion by noting that the data oscillate at the same frequency 
as the first-order interference shown in Eig. 2A and by recalling the discussions of Eqs. (4) and (5). We 
also note that Qi decreases as the degree of the spatial coherence of the source is reduced, making the 
dominant contribution in this case; we reach this conclusion by an examination of Eq. (5), which shows that 
Qii, in contrast to Qn, does not decrease with decreasing degree of spatial coherence of the source. We see this 
behavior in the sequence of results shown in Eigs. 2 (F to H). For example, in Fig. 2F, the contribution from 
C/2^ is smaller than that from Qi. This transition is marked by the formation of second-order correlations in 
the angular position and 0AM variables. 
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Coherent Beam ro = 210 jim ro = 150 jim ro = 70 pm 



Figure 2: Interference transitions in the OAM-mode distribution of light. (A to D) First-order (Young’s) interference. (E 
to H) Second-oder HBT interference. The first column (A and E) shows interference produced by coherent light, whereas the 
other panels show the measured interference for different strengths of the fluctuations of pseudothermal light, as characterized 
by the Fried coherence length. In each case, the angular width of the slits o is 7r/12 and the angular separation of the slits ipo 
is 7r/6. Bars represent data, whereas the line is the theoretical curve predicted by theory. 


It is interesting that there is a regime of random fluctuations for which strong frequency-^ oscillations 
are seen in the first-order interference while strong frequency-2£ oscillations are seen in the second-order 
interference (see Figs. 2, B and F). Note also that, for the case of quantum correlations, entangled photons 
do not produce interference in singles but only in correlations such as those shown in Fig. 2 (D and H) [8p2 


The interplay between Qi and Q 2 i might be useful to the study of the relationship between coherence and 
the quantum nature of light. 


It is important to remark that different degrees of coherence define regimes of the HBT effect 28 


shown in Fig. 2. In our case, the varying relative magnitude of the three terms contributing to the second- 
order coherence results in different shapes (see Eq. 3). For example, makes the pattern in Fig. 2E 
sharper, but the same term changes the frequency of the interference structure in Eig. 2H. 

The general form of the azimuthal HBT effect is obtained when the intensity correlations are calculated 
for arbitrary mode indices ii and £ 2 - As discussed above, the HBT effect depends on the degree of coherence 
of the source. Specifically, an interesting feature is observed for the partially coherent regime characterized 
by tq equal to 150 ^m. In our experimental study of this situation, we hold the 0AM value measured in 
one arm of our interferometer fixed at the value io whereas we vary the 0AM value in the other arm. We 
set the value of £o first to +2 and later to —2. In the other arm, we perform measurements for each value in 
the range £ = —15 to -|-15. The results of these measurements are shown in Eig. 3. It should be noted that 
the 0AM spectrum plotted as a function of the the 0AM value of arm 2 is shifted left (see Fig. 3A) or right 
(Fig. 3B) depending on the value of 0AM chosen for arm 1. The procedure used in the measurement is 
analogous to using one fixed detector and one moving detector in the original setup of HBT . The results 
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of Fig. 3 (A and B) are described by the quantity {hlio) and can be expressed in terms of five contributing 
terms (see section 2 of the Supplementary Materials). 

For the strength of fluctuations that we used for these measurements, one of the detectors measures an 
interference pattern equal to the one shown in Fig. 2C, whereas the other measures a noisy but constant 
signal. When the correlation of the two signals is calculated, the visibility of the interference pattern is 
dramatically increased and shifted in the OAM-mode distribution of the field. Effectively, we are using the 
random fluctuations of the field to increase the fringe visibility. For example, if instead of projecting an 
0AM value equal to 2 or -2 as we did, we could project on i equal to zero and retrieve the original but 
improved pattern with increased visibility. This effect could find importance in realistic applications. These 
effects manifest the presence of second-order correlations in the 0AM components and angular positions of 
pseudothermal light. 


Arm 1 fixed at Iq = —2 Arm 1 fixed at Iq = +2 



Figure 3: Experimental demonstration of the azimuthal HBT effect of light. (A and B) AG^ plotted as a function of the 
0AM value of arm 2 for two different values of the OAM number of arm 1. The green bar shows the center of the interference 
pattern for singles counts shown in Fig. 2C, whereas the purple bar shows the center of the displayed interference pattern. 


We would like to emphasize that although the angular slits and the forked holograms for OAM projections 
are realized on the same SLMs, they correspond to conceptually distinct components of the experiment. The 
angular slits are used to provide a nontrivial azimuthal structure for the incoherent source, whereas the 
forked holograms are used to measure correlations in the OAM domain. 



























































































Measurement of angular momentum correlations and angular po¬ 
sition correlations 


Now we explore the nature of the underlying fluctuation-induced correlations in 0AM and in angular position 
that lie at the origin of the HBT effect. The superposition of randomly fluctuating waves produces an 0AM 
spectrum that broadens with the degree of fluctuation in the source of pseudothermal light. In the present 
experiment, the 0AM spectrum is controlled by setting tq equal to 70 p,m. This situation produces a broad 
0AM spectrum that remains almost constant over the range of 0AM values that we measure. We use the 
same setup as that of Fig. 1, although we omit the two angular slits that we used in the studies of azimuthal 
HBT interference effects reported above. On the first SLM (see Fig. 1), we display a forked hologram 
corresponding to a fixed value of 0AM, whereas on the second SLM, we display a series of holograms with 
different values of 0AM. The measured intensity for a single value of 0AM (7^) that is projected out using 
the SLM can be approximated as / r'^dr'^dcj)^ \£{r)^ gijY, where g{r) is the Gaussian mode supported by 
the SMF (see section 3 of the Supplementary Materials). 

In Fig. 4A, we plot the measured value of /{If,^){Ii^). We find a strong positive correlation 

between the 0AM values measured in the two arms. As shown in section 3 of the Supplementary Materials, 
in the limit of a strong fluctuations, second-order correlations in the 0AM degree of freedom can be described 

by 

{IM = {h,){h,){l +5, ( 6 ) 
Our experimental results show crosstalk between different 0AM numbers that is not predicted by Eq. 
([^. This crosstalk results from experimental imperfections in the projective measurement process used to 
characterize 0AM. The correlations in Fig. 4A show two significant differences from the quantum correlations 
observed in spontaneous parametric down conversion (SPDC). The first is that SPDC shows strong anti¬ 
correlations of the the two 0AM values. This behavior is a consequence of the conservation of 0AM in a 
parametric nonlinear optical process. The second difference is the presence of a background term (the “1” 
in Eq. (®), which prevents the existence of perfect correlations. 

Randomly fluctuating beams also produce correlations in angular position. These correlations are inves¬ 
tigated by encoding angular apertures onto the SLMs. In order to make our measurements precise, we utilize 
narrow angular apertures of 7r/15 radian size. We keep one aperture at a fixed location, and we measure 
correlations for 60 different angular positions of the other aperture. Because of rotational symmetry, this 
procedure permits the full characterization of correlations in angular position. As shown in Fig. 4B, for 


9 


a 


O) 

-j-i 

O 

d) 

'o' 

Sh 

Dh 


for 0AM Correlations 



-10 -5 0 5 10 


Projected £ in Arm 1 

I. 

1.50 1.55 1.60 1.65 1.70 1.75 1.80 1.85 1.90 1.95 2.00 


for Angular Position Correlations 



Figure 4: Measurement of intensity correlations in the angular domain for random light. (A) Normalized second-order correlation 
function in the OAM domain. (B) Presence of strong correlations for the conjugate space described by the angular position 
variable. 


this level of fluctuation, the intensities of the projected angular apertures are strongly correlated, and the 
nature of these correlations can be approximated by + f{(p — ipo)). In this equation, 

the subscript 4> refers to the arm where the variable-position angular aperture is placed and represents 
the arm with the fixed-position aperture. Also, f{4>~ 4’o) represents some strongly peaked function (see the 
Supplementary Meaterials). 

As we have shown throughout this paper, the HBT correlations of pseudothermal light lead to effects 
that show resemblance to those previously observed with entangled photons [^|^[^[^. The reason for 
this behavior is that, in contrast to the degree of second-order coherence that describes coherent light, 
the functions that describe second-order correlations in angular position and OAM for random fields are 
nonseparable. For example, Eq. 6 does not contain the product of the averaged intensities measured by each 
of the two detectors. The presence of a term that describes point-to-point correlations (in this case, the delta 
function 5i^ ) does not allow the factorization of the degree of coherence as the simple product of intensities 

between the two arms. As a consequence, the HBT structures are also described by a nonseparable function, 
and its frequency, visibility, and shifts increase with the fluctuations of the source or the strength of angular 
position and OAM correlations. As the strength of the fluctuations decreases, the nonseparable part of the 
function tends to vanish, and thus, the second order correlation function can be factorized in terms of OAM 
or angular position. A separable function will not lead to the HBT effect in the OAM-mode distribution of 
light; see the transition shown in Fig. 2. 

Intensity correlation in the OAM components and angular position of pseudothermal light show simi- 
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larities with the azimuthal EPR effect, observed in photons entangled in angular position and 0AM [^. 
However, our results show that for pseudothermal light, the correlations are present but not perfect, unlike 
the case of entangled photons where the correlations are perfect. Thus, it is impossible to violate, for ex¬ 
ample, the azimuthal EPR criterion (A£)^(A(/))^ > 1/4. However, as shown in Fig. 4, our correlations are 
stronger for same values of 0AM or angular positions. For example, if background subtraction is performed, 
the variance product for A£ and Ac/ is similar to that achieved for nonclassical light. For our experimental 
results (A€)^(A(/)^ is 0.054, of similar order to the one reported by Leach et al. [^. The uncertainties were 
measured by performing a least squares fit of the data to a Gaussian distribution and recording the standard 
deviation of the result. Note that this does not imply a violation of the EPR criterion. 


Discussion 


The azimuthal HBT effect unveils fundamental physics that can be applied to develop novel applications that 
exploit 0AM correlations in random light. We believe that many interesting protocols for remote sensing and 
object identification that use azimuthal correlations in entangled photons will be able to exploit azimuthal 


correlations in random light and the azimuthal HBT effect 29 -31 . Furthermore, in recent years, researchers 
have developed interest in utilizing beams carrying 0AM for applications in astronomy, but unfortunately the 


propagation through random media produces chaotic phase fluctuations and optical vortices 31-35 . These 


effects pose serious problems for methods based on 0AM of light, limiting their applications . However, 

it has been shown that second-order interference effects are less sensitive to the coherence properties of the 
source. This is one of the advantages of the HBT interferometer against the Michelson interferometer [3^ . 
In addition, it has been demonstrated that imaging schemes based on second-order correlations are robust 


against turbulence 38 . Therefore, we suggest that the azimuthal HBT effect offers the possibility of exploring 


novel phenomena in astrophysics, one example being the relativistic dynamics produced by rotating black 


holes 21 


We have demonstrated that random fluctuations of light give rise to the formation of intensity correlations 
in the 0AM components and angular positions of pseudothermal light. These correlations are manifested 
through a new family of interference structures in the OAM-mode distribution of pseudothermal light that 
can be described by the azimuthal HBT effect. We have shown how the strength of the random fluctuations 
of light determines various regimes for this effect. In addition, we identified two key features of the azimuthal 
HBT effect. The first is characterized by a structure in which the 0AM frequency is doubled with respect 
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to the interference produced by a coherent beam of light. The second is marked by a shift of the 0AM 
spectrum with a change in the 0AM reference value. We anticipate that these properties of random optical 
fields will be fundamentally important for applications where quantum entanglement is not required and 
where correlations in angular position and 0AM suffice. 
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Materials and Methods 


Source of Pseudothermal light. Pseudothermal light was generated by means of phase screen holograms 
obeying Kolmogorov statistics. Kolmogorov’s statistical theory is used to model chaotic turbulent fluids. 
We have generated Kolmogorov phase screens for varying levels of simulated randomness by using the 
approximate power spectral density of d)(/) « 0.023rQ Fried’s parameter rp is related to the 

average coherence length between two points in the beam. By adjusting Fried’s parameter tq, we can 
increase or decrease the size of and the distance between the phase cells and thus the amount of randomness 
in the phase screens. By adding normally distributed deviations to the power spectral density, we can then 
take the real part of the inverse Fourier transform in order to generate a single Kolmogorov phase screen. 

A DMD can be used to manipulate both the phase and amplitude profile of a light beam. A translation 
in a binary diffraction grating will cause a phase shift to occur in the diffracted light, whereas varying the 
duty cycle of the periodic grating will change the efficiency, and thus the amplitude, of the diffracted beam. 
Both of these techniques can be done locally to spatially control the phase and amplitude of the beam. 
The generated Kolmogorov screens were then converted into binary diffraction gratings to be displayed on 
a DMD. 



Figure 5: Example of a frame sent to the DMD. It contains 24 binary holograms encoded in bit plane slices. 


We used the Texas Instruments LightCrafter Evaluation Module (DLPC300) which drives a Texas Instru¬ 
ments DLP3000 DMD. The DMD contains an array of 608 x 684 micromirrors with a total diagonal length 
of 7.62 mm. The DMD was operated in a mode that allowed a binary pattern to be displayed at a rate of 
1440 Hz. The DMD takes a 24-bit color 60 Hz signal over an HDMI connection. Because the image contains 
24 bits, a single video frame can contain 24 binary images. In this mode, the DMD will cycle through the 
least significant bit to the most significant bit in the blue signal of a frame. Then, the DMD will display 
the bits in the red signal and, finally, the green signal. Kolmogorov screens (72,000) were encoded into three 
thousand 24-bit frames for each value of Fried’s parameter tq = 70 /rm, 150 ^m and 210 /xm. Figure 5 
shows an example of one of the generated frames sent to the DMD. This frame contains 24 binary holograms 
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encoded in the bit planes of the image to be displayed sequentially. 



Figure 6: Example DMD holograms and resulting beams measured before the image plane. 


Figure 6 shows examples of the intensity distribution for three random beams generated by this method. 


In addition, in each case, an example of one of the holograms used to produce the beam is shown. 


Note that 


the randomness within the beam increases as the value of Tq decreases. 
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Supplementary Materials 


1. The HBT effect for symmetrically displaced modes (£ and —£). 

Here, we derive the equations utilized in the manuscript. We start by describing HBT interference of 
pseudothermal light. We assume that the beam of light from our laser is described by the electric field £{r), 
where r is the radial coordinate in the transverse plane. In addition, we assume that the initial electric field 
does not possess any azimuthal dependence. We encode a random Kolmogorov phase screen <I>(r, </>) onto the 
beam. Later, the field illuminates two angular apertures centered on angles 0 and (f>o. Thus the field after 
the two slits is given by 


= £{r)<^{r,(j))[A{(j)) + A{(j) - (j)o)].- (7) 

As described earlier in the manuscript, we replace the widely used ground-glass plate with a series of phase 
screens that change rapidly in comparison to the accumulation time of the measurement, thus creating an 
ensemble of field realizations. The next step is to find the intensity of the field for a given 0AM eigenstate. 
We can write the electric field after the two slits as a linear combination of radial-OAM modes. We designate 
a complete radial basis as Rp{r), although we do not make use of any explicit form for this basis. We thereby 
express the field after the slits as 

___ c>ii<P 

( 8 ) 


^ {r, </>) = ^ apiRp (r) 

i,p 


where the expansion coefficients api are given by 

Opf = J rdrdipRp (r) (r, (j)). 

Thus, the measured intensity after projecting the beam onto 0AM mode i is given by 


(9) 


7 ^ = = yi / ridrir2dr2d4>id4>2'^*{ri, (j)i)Rp{ri] 

p p 




f gj^(01-02) 

= / ridrir2dr2d(j)id(j)2'^*{ri,4>i)'^{r2,(l)2) - ^^ Rp{ri)R*{r2) 


( 10 ) 


= Ty rdrd(j}id4>2'i’*{r,(j)i)'ii{r,(j)2)e^^^‘^^~'^^\ 

where the last form comes from using the relation 7?p(r’i)i?*(r2) = (l/ri)5(ri — r2), which is true for any 
complete normalized set of basis functions, where d{r) is the usual Dirac delta function. Now we replace 
'I'(r, (j)) with the electric field after the two angular slits given by Eq. For simplicity, we first approximate 
A{cj)) by 6{(j)) and A{(j) — (j)o) by S{(j) — c/jq). The quantity li then becomes 

^ J rdr|£’(r)|^ {2-f (jio) + ((>o)$(r, 0)} . 


If = 


( 11 ) 
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We next take the ensemble average to obtain 


ih) = ^ y’rdr|£'(r)|"{2 + e-*^^«(<i>*(r,O)<i>(r>o))+e*'^«(<i>*(r,0o)‘f(r,O))}. (12) 

In reality, however, the finite size of the slits produces an envelope, caused by diffraction, that modulates 
the form of the interference pattern. If a slit with a width a is considered, the interference pattern can be 
easily calculated to be , where in our case f{(j)) is equal to 1 in the range from —a/2 to 

a/2 and is equal to 0 otherwise. This integral produces a diffraction envelope given by ^sinc(^). The 
diffraction produced by the second slit can likewise be described as ^sinc(^)e“*^‘^“. The intensity of the 
total diffraction is described as (^)^sinc^(^)(l + cos (t'c/o))■ Taking this result into account, we find that 
the first-order-interference diffraction pattern is given not by Eq. |12| but rather by 


{h) = 


a^sinc^(£a/2) 

27r^ 


J rdr£^{r) {2 -F ($* (r, 0)$(r, </o)) + ($* (r, /)o)$(r, 0))} . (13) 


We next develop appropriate approximations for the quantities defined above. A reasonable assumption 
is that the field fluctuations follow Gaussian statistics such that 

rj -I- r| — 2rir2 cos (j)o . 


($*(ri, 0 )$(r 2 ,(/)o)) = exp(-- 
By setting ri equal to r 2 , we find that 


-)■ 


4^21 sin ^a|2 

($*(r,0)$(r,())o)) = exp(- ^ ^ ) = exp(-/3r2). 


(14) 


(15) 


The last form of this expression defines the quantity /3. For a fully coherent beam, (that is, for rg ^ r, where 
cq is the Fried parameter introduced in Section 3) we see that to very high accuracy (<f>*(r,0)$(r,c/o)) is 
equal to 1. As tq decreases the value of the correlation function ($*(r, 0)$(r, c/o)) also decreases. Through 


use of Eq. 14 expression (6) for the intensity can be expressed as 


{le) =Xq-\-Xi 

Xq ^ ^ J rdrS'^ir) (16) 

Xi = J rdr\£{r)fexp{-l3r^). 

As To decreases, /3 increases and most of the contribution to the integral in the last expression comes from 
r « 0, but since the integrand is zero at that point, the integral vanishes. In this limit, Xq makes the only 
contribution to (/^) and the spectrum becomes flat, that is, (Ii) shows no dependence on the value of i. 
Next we derive an expression for the correlations between projections onto 0AM values of i and that 
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IS, 


^ J ridri|£:(ri)|'{2 + e-*^'^“$*(ri,0)$(ri,</.o) + e*^‘^'>$*(ri,(Ao)$(ri,0)} 

X y'r2dr2|£(r2)n2 + e*^^‘>$*(r2,0)$(r2,<^o) + e-*^^«$*(r2,<^o)$(r2,0)} (17) 

= Qo + Qt + Qii, 

where the three contributions to lil-i are given by 

J ridrir2dr2\S{ri)\'^ \£{r2)\'^ + ^ J ridrir2(ir2 |£(ri)|^ |£(r2)|^ $*(ri, 0)$(ri, 0o)^’*(?'2,0)$(r2, (/>o) 

+ ^ |£(ri)|^ |£:(r2)|^$*(ri,())o)$(ri,0)$*(r2,<?io)^’(?'2,0) 

(18) 


gi=— J ridrir 2 dr 2 \£{ri)f \£{r 2 )\^ e {$*(ri, 0)$(ri, i^o) + $*(?'2, ()'o)^’(?'2,0)} + c.c., and (19) 

G 2 e J ridrir 2 dr 2 |f(ri)|^ |£’(r2)|^ {$*(ri, 0)$(ri, ())o)$*(r2, ^o)®(?’2,0)} + c.c. (20) 

We next estimate the ensemble averages of these quantities. For a field with strong random fluctuations, 
the field correlation between two different angular positions is very small and thus the term G^ does not 
contribute significantly to A similar situation occurs for the second and third contributions of Go] 

it is important to note that these terms contain the quantities ($*(ri, 0)$(ri, ())o)$*(r2, 0)<i)(r2, (?iio)) and 
(<i)*(ri, ^o)‘I’(?'i) 0)'h*(r2, ^io)‘J^(’'2,0)), and these quantities vanish when ri = r 2 . They vanish because they 
describe the average of the product of two chaotic and independent variables. Thus, the main contributions to 


the second-order interference are the first term in Eq. 18 (which does not vary with (j)Q) and the contribution 

„ —2i^0o 


{G2e) =- 


dTT^ 


ridrir 2 dr 2 \£{ri)\^ \£{r 2 )\^ ($*(ri, 0)$(ri, ^o)®*(?'2, ()'o)^(?’2, 0)) +c.c.. 


( 21 ) 


It is important to note that, contrary to the correlation functions given by (<F*'(ri, 0)$(ri, ())o)<i)*(r2, 0)$(r2, (/>o)) 
and ($*(ri, (/)o)$(r'i, 0)<i)*(r2, 0o)‘l’(?'2,0)), the correlation function ($*(ri, 0)‘l>(ri, ^o)‘l’*(^2, </'o)‘&(^2, 0)) is 
equal to unity for ri = r2. We therefore obtain 


{G 21 ) = 


47r^ 


r^ciri \£{ri)y + c.c. 


( 22 ) 


We thus conclude that the quantity (IiI-i) is given by the sum of the contributions of Eq. 22 and the first 
term in Eq. [T^ or 


(///_;)« J rdr|£:(r)|^^ -b J r‘^dr\£{r)\‘ 


(23) 


For the case of slits with finite size, this result must be modified for the same reasons given in the discussion 
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following Eq. 12 One thereby obtains 


(/,/_,)« Jrdr\£{r)f^ + Jr^dr\£{r)\‘ 


(24) 


2. The HBT effect for arbitrary mode indices ii and £ 2 - 

The intensity correlation between two arbitrary 0AM modes produces a complicated second-order 

correlation function comprised of five terms. The contribution of each term is determined by the degree of 
fluctuations in the field. One is the constant term Oo given by / ridrir 2 dr 2 {\£{ri)^ |£’(r 2 )|^). There are 
two terms whose contributions are equally important; one oscillates with a frequency £i4>o and the other 
with £ 2 (t>o- The strength of these terms is determined by the quantity ($*(r, 0)$(r, which is negligible 
for highly chaotic light (rg r). The frequency of the fourth component is determined by the quantity 
(£1 -I- £ 2 ) 00 , although its strength is dictated by the quantity (<i)*(r, 0)<l>(r, ^ 0 )^ 1 ?* (»') 0)d>(r, ^o))- For highly 
chaotic fields this is an extremely small contribution. The primary contribution to ) is therefore given 


by the term 
Qild2 = 


4 • 4 

a Sine 


^ [^ 1 — 


47r^ 


ridrir 2 dr 2 |f(ri)|^ 1^(?' 2 )|^ (e 0)$(ri, (/)o)4>*(?’2, (()o)$(?' 2 ,0)) 


-l-c.c.) 

(25) 

Note that this contribution describes an interference pattern that depends on the values of both £1 and £ 2 - 
3. Orbital angular momentum correlations and angular position correlations. 

In this section we derive expressions for the correlations of pairs of 0AM values and pairs of angular 
positions. The light that emerges from the DMD is given by £(r)$(r, </)). We make two copies of this field 
using a beam splitter and find the coincidences between projections onto two different modes of light. Let 
us first discuss the projection of one of the beams. The amplitude of the projection onto 0AM mode £ is 
given by 


ai = y rdrd())f(r)4>(r»--^g(r), 


(26) 


where g{r) is the radial profile of the single-mode collection fiber, which is a Gaussian function. The intensity 
is given by 

n pit(p2 

ridrid(l)i£{ri)^{ri,(j)i)—j^g{ri)x J r2dr2d(j)2£*(r2)‘^*{r2,(t)2)^^j^g{r2). (27) 


Therefore the ensemble-averaged intensity after the projection is given by 


(le) = / ridrir2dr2d(j)id(j)2£(ri)g(ri)£*(r2)g(r2)- 


o-i^4>i-4>2) 


27r 


-($(ri,(()i)$*(r2,())2))- 


(28) 


We are considering the case of highly fluctuating light; in this regime ($(ri, ())i)<i)*(r' 2 , ^ 2 )) can be approxi- 
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mated by {\/ri)5{r2 — ri,(j )2 — leading to the result 

^ J ridridf/ji \£{ri)f g{riy. 


(29) 


Note that this quantity is independent of the value £. Now let us consider the case of two coincident 

projections. The amplitude of coincident projections is given by 

/ r g-iti'Pi \ 

= / r,drid(j)i£{ri)^{ri,(j)i) — j^g{r,) . 


We measure the rate at which these coincidences occur, which is given by 


leje^ = ridrid(j)i£{ri)^{ri,(l>i)- 


Mi4’3 


- 5 (^ 1 ) X / r2dr2d4>2£{r2)^(r2,(j)2)- 


if 2 04 


■ 9 (^ 2 ) 


rzdr'id(j)^£*{r^)^*{r'i,(l)3,)—j=g{r'i)x J r4dr4d(l)4£*{r4)^*{ri,(l)^)-j=g{r4). 


(30) 


(31) 


After taking the statistical average, we obtain 


{hjt2) = / ridrid4>i£{ri)- 


,if 103 


■9{ri) X / r2dr2d(j)2£{r2)- 


o — i^24>2 

,2^2 <p4 


- 5 (^ 2 ) 


/ g 2 £l 03 r g*t 2 94 

r 3 dr 3 d(j) 3 £*{r 3 )-j=g{r 3 ) x J ridrid4>i£*(ri)-j^g{ri) x (4'(ri, ^i)$*(r3, ())3)4>(r2, ^2)$*(f4, (('4))- 


(32) 


Following the same considerations and a similar procedure those used in obtaining Eq. 17 we write the 
four-point coherence function as the sum of three contributions, each a product of two 2-point coherence 
functions. One of the contributions is always negligible for highly chaotic light. Another contribution leads 
to the simple product {Ig^){Ii 2 )- This contribution is actually independent of the values ii and £2 in the limit 
of highly chaotic light for the same reason stated above in connection with Eq. The last contribution is 
given by 


ridrid(j)i£{ri)- 


„-ifi0i P gif204 

g(fi) X / r4dr4d(j)4£*{r4) g{r4) x {^{ri,(j)i)^*{r4,4>4)) 


.—i£24*2 


pi^l<p3 


(33) 


/C — It2<p2 r 

r2dr2d4>2£{r2)—j^g{r2) x J r3dr3d4>3£*(r3)—j=g{r3) x (4'(r2, ^2)$*(f3,'('a))- 
By invoking the same approximation used above to evaluate the coherence functions as delta functions, we 


simplify this expression to 


rdrd(j) \£{r)^ di^Y 


o — 2 (^ 1 —^ 2 )^^ 


27r 


(34) 


Note that the integral over </> vanishes unless £i = £ 2 - Thus, we finally obtain 


{^ 1 / 2 ) = (-^fi)(-ff2)(l +'^fi.f2). (35) 

which is the expression used in the body of this paper to explain our experimental results. Note that the 
correlations between two different values of 0AM are half as large as those between the same value of 0AM. 
We can perform a similar calculation to find the correlations between two angular positions. For the case 
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of a single beam (no beam splitter), the amplitude of the projection for a single value of (j) is given by 

00 = y rdrE(r)^(r,(j))g{r). (36) 

It follows that the ensemble-averaged intensity at one of the detectors is given as follows 

{l4>) = M (37) 

= ( J ridri£{ri)^{r,I (j))g{ri) x J r 2 dr 2 £*{r 2 )^*ir 2 , (I))g{r 2 )) (38) 

= J ridrir2dr2£{ri)£*{r2)giri)g{r2) {^{r,I (I})^{r2,(lj)*) (39) 

= J rdr\£{r)f g{r)'^. (40) 

The last form follows from the substitution (<I’(r,i (f>)^{r 2 , (j))*) = (l/ri)(5(ri — r 2 ). If we now add the beam 
splitter and find the probability for coincidence detection of two beams, we obtain 

{U 1 U 2 ) = J ridri£{ri)g{ri) X J r2dr2£{r2)g{r2) (41) 

X J r 3 dr 3 £*{r 3 )g{r 3 ) X J r4dr4£*{r4)g{r4) x {^{ri,(l>i)^*{r3,(l)3)^{r2,(l>2)^*{r4,(l>4))- (42) 

Similar to the 0AM case, we find that this expression yields two contributions. One of these contributions 
is independent of the values (pi and 4>2 whereas the other contribution is large only if the two intensities are 
measured at pi = (j) 2 - This result can be described by the relation 

{IM = (/ 0 ,)(/ 0 ,)[l + 5{4>4 - cP2)]. (43) 

As expected, and similar to the case of correlating two OAMs, two non-overlapping angles share no corre¬ 
lation. However, in contrast to the 0AM variable, the angular position variable is not discrete, and one is 
allowed to correlate two regions defined by two angular positions that are not orthogonal, and consequently 
there is a partial overlap between the two correlated regions. Thus the degree of second-order correlation 
can take any value between 1 and 2. Therefore an appropriate expression for this correlation function is 
given by 


(701— (/0i)(702)[l + f{4>i - ((> 2 )], 

where /(</>! — (^ 2 ) is defined as 

{ 0 for 1^1 — (p 2 \ > a /2 

^_| 0 i^-^ for \(j)i - (p 2 \ < a/2, 

/{(pi — P 2 ) can be interpreted as the fractional angular overlaps of the two slits. 


(44) 
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